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Abstract. In recent years, it has become clear from multidimensional simulations that 
the outcome of deflagrations depends strongly on the initial configuration of the flame. 
We have studied under which conditions this configuration could consist of a number 
of scattered, isolated, hot bubbles. Afterwards, we have calculated the evolution of 
deflagrations starting from different numbers of bubbles. We have found that starting 
from 30 bubbles a mild explosion is produced (M( 56 Ni) = 0.56 Mq), while starting 
from 10 bubbles the star becomes only marginally unbound (K = 0.05 foes). 


1 The initial configuration of deflagrations in white 
dwarfs and its outcomes 

The explosion mechanism(s) responsible for thermonuclear supernovae (SNIa) 
is still not well known. Although there have been recently some claims that the 
delayed detonation mechanism lacks a physical background, there are still un¬ 
explored mechanisms by which the transition from deflagration to detonation 
could occur (see Garcia-Senz & Bravo, this volume). Here we concentrate on 
the other possibility, i.e. a pure deflagration that would process about a solar 
mass, synthesizing of the order of O.5M 0 of 56 Ni in order to make a typical 
SNIa. Multidimensional calculations of deflagrations bear the advantage over 
one-dimensional models that the energy generation rate becomes eventually in¬ 
dependent from the local value of the flame speed (|^|0, see also Niemeyer, this 
volume). In these calculations, the acceleration of the fuel consumption rate is 
due to the deformation of the flame surface, which is well accounted for by 3D 
hydrocodes. 

In this picture, the main free parameter is the initial configuration of the 
flame. There have been a few || simulations of the transition from the hydro¬ 
static phase to the hydrodynamic one, but a self consistent multidimensional 
initial model is still far from being available. However, a plausible ignition sce¬ 
nario was suggested in || : the nearly-simultaneous runaway at several different 
spots (from here on, bubbles) in the central region of the white dwarf. In this 
paper, we explore the dependence of the outcome of the explosion on the initial 
number of bubbles. First, we address the statistics of the initial distribution of 
bubble radii and, then, we present the results of a couple of hydrodynamical 
simulations performed with an SPH code (@ and Garcfa-Senz & Bravo, this 
volume). 
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1.1 Statistical approach to the initial distribution of igniting 
bubbles 

Our statistical approach is based on the following assumptions: A) As a result 
of convection, there appear a number of hot spots, characterized by its central 
(peak) temperature, To, and its thermal profile (which we take here exponential, 
with characteristic exponent Rq )• B) Each hot spot evolves adiabatically in 
place, in a time given by its ignition timescale. C) The peak temperatures of 
the hot spots can be characterized by a statistical continuous function, 9 (To). 

Depending on the thermal gradient inside each hot spot our results can be 
split into two different regimes: 1) If the thermal profile is shallow enough, 
bubbles grow due to spontaneous flame propagation, 2) otherwise, bubbles grow 
conductively. In the first case, the distribution function of bubble radii is a time 
function given by 

dA _ A 1 ' 3 ■ 10 9 Tj (T ob ) exp (BR h /R 0 ) 

di? b ° Ro {i + Ti (To;,) [1 — exp (SI? b /i?o)]} 1+1/B 

where t; (T 0 b) ~ A (T 0 f,/10 9 K) B is the ignition timescale, and T 0 b = 10 9 K. In 
the second case, the distribution function becomes 


In the first case, and for any reasonable choice of the function 9 , the resulting 
distribution has a sharp peak at a determined value of R b so that the initial con¬ 
figuration is composed of an arbitrary number of equal size bubbles. In contrast, 
in the second case the bubble radii distribution is continuous up to a maximum 
radius, which results in a initial distribution of unequal size bubbles. 

1.2 SPH simulations of deflagrations triggered by floating bubbles 

Given the constraints that the value of Rq must satisfy in order to get an initial 
configuration composed of identical bubbles, we estimate as quite improbable 
this case. However, for simplicity, here we have adopted this configuration as the 
starting point of our simulations. We have followed the evolution of the explosion 
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Table 1. Model parameters and results 


Model 

N b 

Rb 

Mb 

in 

Mni 

Mc+o 

Msi 



(l0 6 cm) 


(10 51 erg) 

(Me) 

(M 0 ) 

(M 0 ) 

B30U 

30 

6.7 

2.8% 

0.44 

0.56 

0.65 

0.05 

B10U 

10 

9.7 

3% 

0.05 

0.24 

0.88 

0.04 
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Fig. 1 . Snapshots of the temperature distribution in a meridian plane at times 0, 0.27, 
0.58, and 0.81 s for model B30U. The temperature scale is shown at the bottom left of 
the image, while the length scale is shown at the top left of each snapshot (the length 
of the vertical bar is equivalent to 200 km) 


of a Chandrasekhar mass white dwarf of initial density 1.9 10 9 g/cm 3 starting 
from two different numbers of equal size bubbles, as detailed in Table [I] (iVb 
is the initial number of bubbles, i?b its initial radius, Mb the mass incinerated 
initially, and the other symbols have their usual meanings). The calculation was 
performed in 3D with the above mentioned SPH code using 250,000 particles, 
with a central resolution of 20 km. It is important to emphasize that we did not 
impose any artificial symmetry conditions and, thus, our model has no artificial 
characteristic length other than the own resolution of the code. Actually, the im¬ 
position of artificial symmetry conditions (for instance, symmetry planes) in any 
3D hydrodynamic calculation biases the development of large-scale structures or 
even of the small-scale ones if they are close enough to a symmetry plane. 

The results of our calculations can be found in Table |l| and in Figs. |l] and §a. 
The outcome depends strongly on the initial number of bubbles. The deflagration 
is powered by the evolution and interaction of the bubbles: growth, buoyancy 
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Fig. 2. Results of the hydrodynamical calculation, (a) Nuclear energy generation rate 
as a function of time for models B30U (solid line) and B10U (dotted line), (b) Final 
distribution of elements in velocity space for model B30U 



Fig. 3. Spatial distribution of the main chemical species at the end of the calculation 
of model B30U, in the same meridian plane as in Fig. [IJ 


(second snapshot in Fig. |l|), and merging (third snapshot). The maximum ac¬ 
celeration of combustion is obtained when the bubbles interact with each other, 
feeding a rich spectrum of scalelengths to the hydrodynamic instabilities (third 
snapshot in Fig. [j], see also Fig. [in). It is not a surprise that this interaction is 
favoured in the presence of a large number of bubbles. The distribution of nuclei 
in velocity and space for model B30U is shown in Figs. |]b and [|. Our results 
agree qualitatively with those obtained in || . 

In summary, our best model (B30U) fails to give the magnitudes adequate 
for a typical SNIa. A possible cause is a poor representation of the subsonic 
flame at low densities (a common problem in most SNIa hydrocodes). However, 
a thorough exploration of the parameter space of initial conditions is in order. 
This work has been supported by the MCYT grants EPS98-1348 and AYA2000-1785 
and by the DGES grant PB98-1183-C03-02. 
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